Nanotube field of Ceo molecules in carbon nanotubes: atomistic versus continuous 

tube approach 
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We calculate the van der Waals energy of a Ceo molecule when it is encapsulated in a single- walled 
carbon nanotube with discrete atomistic structure. Orientational degrees of freedom and longitu- 
dinal displacements of the molecule are taken into account, and several achiral and chiral carbon 
nanotubes are considered. A comparison with earlier work where the tube was approximated by a 
continuous cylindrical distribution of carbon atoms is made. We find that such an approximation is 
valid for high and intermediate tube radii; for low tube radii, minor chirality effects come into play. 
Three molecular orientational regimes are found when varying the nanotube radius. 



I. INTRODUCTION 

The discovery of carbon nanotubes (CNTs) by lijima 
Pi and their subsequent large-scale production was 
followed by the synthesis of CNTs filled with atoms 
and/or molecules. These novel hybrid materials often 
exhibit one-dimensional characteristics and are presently 
the subject of fundamental studies as well as research 
aiming at their application in nanotechnology. For a re- 
view on CNTs and their filling we refer to Refs. 3, 4 and 
0, 1^ respectively. Self-assembled chains of Cgo fuUerene 
molecules inside single-walled carbon nanotubes (SWC- 
NTs) , the so-called peapods 7} , provide a unique exam- 
ple of such nanoscopic compound materials, and feature 
unusual electronic [g and structural properties. High- 
resolution transmission electron microscopy observations 
on CNTs filled sparsely with Ceo molecules Q demon- 
strate the motion of the fuUerene molecules along the 
tube axis and imply that the interaction between Cgo 
molecules and the surrounding nanotube wall is due to 
weak van der Waals forces and not to chemical bonds. 

Recently, the way the Ceo molecules of a 
(Ceo)w@SWCNT peapod 0, [13 — Ceo molecules 
inside in a SWCNT — are packed in the encapsulating 
tube has been investigated both experimentally and 
theoretically [H El 111 El IH US- Obviously, the 
structure of a peapod is governed by the interactions be- 
tween the Ceo molecules, and by the way a Ceo molecule 
interacts with the surrounding tube wall. Already 
when considering the stacking of cylindrically confined 
hard spheres, a possible rudimentary description of 
a (C6o)Ar@SWCNT peapod, various chiral structures 
of the s phe res stacking for varying tube radius are 
obtained [Ij • In Ref. E2i Hodak and Girifalco calculated 
lowest-energy (Ceo)Ar@SWCNT peapod configurations 
by means of a continuum approach for the Ceo-tube 
interaction: both a SWCNT and a Ceo molecule are ap- 
proximated as a homogeneous surface — cylindrical and 
spherical, respectively. Although in doing so any effect 
of tube chirality and/or molecular orientation can not be 
accounted for, such a model provides useful information 
about the spatial arrangement of the spherical molecules 
in the tube. Ten different stacking arrangements were 



obtained for the tube radius Rt ranging from 6.27 A 
to 13.57 A. The simplest configuration (Ceo "spheres" 
aligned linearly along the tube axis) occurs for the 
smallest tubes (6.27 A < i^x < 7.25 A). Other phases 
consist of zig-zag patterns or Ceo balls forming helices. 
Some of the predicted phases have been observed exper- 
imentally [iSj. Interestingly, experimental observations 
of similar structures formed by Ceo molecules inside 
BN nanotubes have been reported as well j^. An 
atomistic molecular dynamics study on the arranging 
of Ceo molecules inside SWCNTs was carried out by 
Troche et al. E3| ! the Ceo-tube interaction was modelled 
by adding carbon-carbon Lennard- Jones 6-12 potentials. 
Troche et al. E3 concluded that the chirality of the 
encapsulating SWCNT has only a minor effect on the 
lowest-energy configuration of the Ceo molecules and 
their obtained arrangements, thus depending on the 
tube radius only, are in full agreement with those of 
Hodak and Girifalco [l3- Conclusions on the individual 
orientations of Ceo molecules inside a SWCNT were not 
given by Troche et al. [H — their goal was to study 
the packing of several molecules. Molecular orientation 
effects are expected to come into play at sufficiently low 
temperatures when orientational motion is frozen, and 
indeed do so as was shown in Refs. and 0, where 
the potential energy of a single Ceo molecule confined 
to the tube axis of a SWCNT, called "nanotube field", 
was calculated by treating the tube as a homogeneous 
cylindrical carbonic surface density but retaining the 
icosahedral features of a Ceo molecule. A specific 
dependence on the tube radius was found; three distinct 
molecular orientations were observed within the range 
6.5 < -Rt < 8.5 A. It is our opinion that, for calculating 
tube-Ceo interactions, taking the detailed molecular 
structure of a Ceo molecule into account has priority over 
the chiral structure of a nanotube. Replacing a SWCNT 
by a continuous cylindrical distribution of carbon atoms 
is intuitively justifiable, but treating a Ceo molecule as 
a sphere (as in Ref. [I2j ) with no further structure is a 
more questionable approximation. Indeed, whereas the 
carbon-carbon bonds in a CNT are of one type, a Ceo 
molecule features longer ("single") and shorter ("dou- 
ble" ) bonds, arranged in pentagons — electron-poor 
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regions — and hexagons — electron-rich regions. The 
importance of taking the detailed molecular structure 
properly into account follows from Refs. and 0; but 
the neglect of the discrete atomistic structure of the 
tube when considering Cgo-tube interactions, although 
intuitively plausible, requires solid grounds. The goal 
of this paper is to answer the question how good a 
smooth-tube approximation really is, and to confirm the 
relevance of the precise structure of a Ceo molecule, i.e. 
the importance of allowing for molecular orientational 
degrees of freedom. 

The content of the paper is as follows. In Sec. (nj we 
discuss formulas for the calculation of the nanotube field 
of an encapsulated Cgo molecule for both a "continuous" 
and a "discrete" tube. Then (Sec. 1111(1 . we plot nan- 
otube fields for a selection of representative nanotubes 
and make preliminary visual comparisons between the 
two approaches. In Sec. IIVI we present an all-variable 
treatment and apply it for tubes with intermediate and 
small tube radii. Finally, general conclusions are given 
(SecEJ. 




FIG. 1: A single Ceo molecule (radius -Rcgo) in a SWCNT 
(radius -Rt). Shown is a projection onto the {x, 2)-plane. The 
center of mass of the Ceo molecule is chosen as the coordinate 
system's origin. The tube's long axis coincides with the z- 
ajcis; the Ceo molecule is put in the standard orientation. 



II. NANOTUBE FIELD 

We consider a Cgo molecule in a SWCNT, the molecule 
assuming a centered position in the tube, and set up 
a cartesian system of axes {x, y, z) so that the z-axis 
coincides with the tube's long axis and contains the 
molecule's center of mass (Fig. QJ. The potential energy 
V of the Cgo molecule then depends on the orientation 
of the molecule, which can be characterized by three Eu- 
ler angles (a, /?, 7), on the position of the molecule along 
the tube, i.e. the z-coordinate of the molecular center of 
mass for which we write and on the tube indices 
(n, m): 

y = y(a,/3,7;C;n,TO). (2.1) 

For the Euler angles we use the convention of Ref. 
[1^: a coordinate function /(r = (x, is trans- 

formed as 9^(a, /?, 7)/(r ) = /(9^~^(q!,/3, 7)r), where 
$n(a,/3,7) = 9^^(7)^ny(/3)fR^(a) stands for the succes- 
sion of a rotation over < a < 27r about the z-axis, 
a rotation over < < tt about the j/-axis, and a ro- 
tation over < 7 < 27r about the z-axis again. The 
y- and z-axes are kept fixed. Note that the coordi- 
nate transform associated with the Euler angles reads 
f' = 5H-Ha,/3,7)r = 9l^(-a)fny(-/3)?l2(-7)r and that 
the rotation of the Cgo molecule over —a about the 
z-axis is performed last. As the starting orientation 
[(a = 0, = 0, 7 = 0)] wc take the so-called standard 
orientation [Fig.|2Ia)]: twofold molecular symmetry axes 
then coincide with the cartesian axes and every carte- 
sian axis intersects two opposing double bonds. (We re- 
call that the carbon-carbon bonds of a Cgo molecule can 
be divided into two categories: 60 single bonds, fusing 
pentagons and hexagons, and 30 double bonds, fusing 



hexagons. The latter are somewhat longer than the for- 
mer [13.) Bearing in mind the results of Refs. Hsl and fl^ 

and anticipating the results obtained in the present work, 
we point out two more molecular orientations of impor- 
tance. The first is the "pentagonal" orientation class, 
obtained by the Euler transformation (a arbitrary, fi — 
cos~^ —, — 2 K, 58°, 7 = 0), resulting in two opposing 

VlO+2%/5 

pentagons of the Cgo molecule being perpendicular to the 
z-axis [Fig. |2Ib)]. The second is the category of "hexag- 
onal" orientations, a result of the Euler transformation 
(a arbitrary,/? = cos^^ ^ 21°, 7 = 0), making two 

opposing hexagons lie perpendicular to the z-axis [Fig. 
M,c)]- The angle ^0 = cos"^ is related to the dihe- 

dral angle ip (the inner angle between adjacent faces) of 
a regular icosahedron: ip — tt — 2/3o. Other (/3,7) pairs 
yield "pentagonal", "hexagonal" and "double-bond" ori- 
entations as well: 12 pairs correspond to a "pentagonal", 
20 pairs to a "hexagonal", and 30 pairs to a "double- 
bond" orientation since a Ceo molecule has 12 pentagons, 
20 hexagons and 30 double bonds. 

For the description of the interaction between the Cen 
molecule and the nanotube we follow earlier work |20| 
and treat the Ceo molecule as a rigid cluster of interaction 
centers (ICs). Not only C atoms ('a') act as ICs, but also 
double bonds ('db') and single bonds ('sb'). We label the 
60 atoms by the index Aa — 1, . . . ,60. In the center of 
every of the 60 single bonds an IC is put, labelled by 
the index Asb = 1,...,60. On each of the 30 double 
bonds, 3 ICs dividing the bond in four equal parts are 
put, totalling to 90 db ICs, labelled Adb = 1, . . . , 90. Such 
a construction was originally introduced for modelling 
intermolecular interactions in solid Ceo (Ceo fuUerite); 
having three ICs per double bond reflects the electronic 
density being smeared out along a double bond . 
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FIG. 2: Special orientations of the Ceo molecule: (a) standard orientation (a = /3 = 7 = 0, double bonds perpendicular 
to the 2-axis), (b) "pentagonal" orientation (pentagons perpendicular to the z-axis), (c) "hexagonal" orientations (hexagons 
perpendicular to the z-axis). Double bonds are drawn bolder than single bonds; dashed lines are located "beneath" the {x,y)- 
plane (2 < 0). The surrounding tube is represented as a circular projection onto the {x,y) -plane. Note that a rotation about 
the z-axis over —a does not change double bonds, pentagons or hexagons — for Figs, (a), (b) and (c), respectively — being 
perpendicular to the z-axis. 



Every IC of the Cgo molecule interacts with every atom 
of the nanotube via a pair interaction potential w'(|rv — 
taJ), depending on the type of IC ('t' — 'a', 'db', 'sb'). 
The total potential energy is then obtained by summing 
over all pair interactions: 

^ = E E Y.'^W^r-rxA): (2.2a) 

T 't'— 'a'.'db'.'sb' At 

where r indexes the atoms of the tube and fv stands for 
their respective coordinates. As in Refs. and E we 
use Born-Mayer-van der Waals pair interaction poten- 
tials: 

t,t(r) = ^. (2.2b) 

Again, the use of such pair potentials was originally intro- 
duced for studying Ceo - Ceo interactions in Ceo fuUerite 
|20| : it lead to a crystal field potential and a structural 
phase transition temperature |2lL |2^ in good agreement 
with experiments. The pot ential constants C*, C2 and 
used are those of Ref. Hm. In Eq. H2.2a|l . the sum over 
tube atoms, labelled by the index r and having coordi- 
nates TV = {xr,yT, Zt), can be restricted to atoms in a 
certain vicinity of the Ceo molecule, realized by imposing 
the criterion 

-^min — Zf ^ ^max; (2.2c) 

with Zmin and z^ax cut-off values ensuring convergence. 

In Refs. 10 and 0|, a smooth-tube approximation to 
Eq. (|2.2all was presented. The actual network of carbon 
atoms making up the SWCNT is replaced by a homoge- 
neous, cylindrical "carbonic" surface density with value 
a (units A~^). The Ceo molecule-nanotube interaction 



energy is then rewritten as 

Smooth 

^aRT dZ E«'(Ip-^^aJ), 

't'='a','db','sb' At 

(2.3) 

where p ~ {Rt,^,Z) is the cylindrical coordinate of a 
point on the tube {x — i?TCOs$, y — i?Tsin<I>, z = Z) 
and i?T is the tube radius. The motivation for intro- 
ducing approximation (|2.3|) is twofold. One reason is the 
dependence of Kmooth on the tube radius i?T rather than 
on the tube indices (n,m). Indeed, Rt remains the only 
relevant tube-characteristic parameter and as such sim- 
plifies a systematic investigation of carbon nanotubes. A 
further consequence of the tube's cylindrical symmetry is 
the irrelevance of the Euler angle a (a final rotation of the 
Ceo molecule over —a about the tube axis doesn't mat- 
ter) and of the z-coordinate C (for infinite or long-enough 
tubes). A second advantage of the smooth-tube ansatz 
is the possibility of performing an expansion of V^mooth 
into symmetry- adapted rotator functions, a point we will 
return to in Sec. We stress the limited dependence of 
V^smooth by writing 

^smooth = V;niooth(^,7;-RT)- (2.4) 

To distinguish the smooth-tube approximation from the 
discrete case, we add the subscript 'discrete': 

V = Vdiscrctcia,l3,j;C;n,m), (2.5) 

where the actual expression is given by Eqs. - 

In this paper we test the validity of smooth-tube ap- 
proximation (|2.3|l by comparing Vdiscrctc and Vsmooth for 
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a selection of tubes. Bearing in mind the three qualita- 
tively different radii ranges {i?T < 7 A, 7 A < i?T < 7.9 
A and 7.9 A < Rt) obtained in Ref. 0|, we have selected 
zig-zag, armchair and chiral tubes with radii around 
Rt = 6.5 A, Rt = 7.5 A and Rt = 8.5 A. We have gen- 
erated (n, to) tubes starting from a graphene sheet with 
basis vectors ai = aex and 02 — a^ex+a^ey, where ex 
and ey are planar cartesian basis vectors, and performing 
the roll-up along the vector C(n, to) = nai +md2 [3.l23j|. 
The tube is then positioned so that the C atom origi- 
nally (before rolling up) at Oex + OeV lies in the (x, y) 
plane with x-coordinate and y-coordinate Rt and that 
the cylinder containing the C atoms has its long axis co- 
inciding with the z-axis. The Cgo molecule is initially 
positioned so that its center of mass lies at the origin 
(C = 0); a translation along the z-axis away from the 
initial position is measured via the center of mass' z- 
coordinate C. The radius of the tube with indices (n, to) 
reads Rt = ^Vn'^ +nm + m'^, with a = 2.49 A ^E^: 
the corresponding surface density has the value 



^/3a2 



0.372 A-^ 



(2.6) 



A further tube parameter is its translational periodic- 
ity Az, relevant when considering the ^-dependence of 
Vdiscrctc- While Az is small for non-chiral — i.e. zig- 
zag, Az = VSa, and armchair, Az = a — tubes, the 
translational period can get very large for chiral tubes 
0. A tube may also have an s-fold symmetry axis (co- 
inciding with the z-axis) and therefore a rotational pe- 
riod Aa = 27r/s. When considering a tube with s-fold 
rotational symmetry it suffices to examine the interval 
< a < Aa. The periodicities and other tube charac- 
teristics of our selected tubes are listed in Tabled 



III. MERCATOR MAPS 

To get a preliminary idea of how T^mooth and Vdiscrctc 
compare, we have simply plot t4mooth(/3, 7, i?T) and 
Vdiscrctc(a = 0,/3, 7;C = 0;n,m) for each of the selected 
(n, to) tubes in the form of Mercator maps |2J| . We stress 
that (a = 0, C = 0) is but a particular case and that final 
conclusions should be made not only on the variation of f3 
and 7 but on the varying of a and as well, as we will do 
later on. We do point out, however, that we expect the 
a- and C-dependencies to be of a lesser magnitude than 
the f3- and 7-dependencies since the former correspond to 
a (final) rotation of the molecule about the z-axis (tube 
axis) over —a and a translation of the molecule along the 
z-axis, respectively, and hence relate to the tube struc- 
ture rather than to the molecule structure. (As argued 
in the Introduction, a carbon nanotube can be regarded 
as being more "continuous" than a Cgo molecule.) 

As for the cut-off values, we have found — for C = 
— that Zmin = — 50 A and z^ax = 50 A yield sufficient 
convergence. Note that the choice of Zmin and Zmax fixes 
the numbers Nt of atoms to be taken into account in 



sum (|2.2a|l . In principle, the tube fragment of length 
L — Zmax — ^min has a surface density a = Nt/ {2ttRtL), 
differing from a. We observe that differences between a 
and a are small, however. Although possibly (slightly) 
enhancing the agreement between Vdiscrctc and V^niooth, 
we have chosen not to calculate Vsmooth with a since 
it somehow relates to the tube structure — Nt de- 
pends on (n, to) — hence surpassing the smooth-tube 
approach's underlying basic idea (i?T-dependence rather 
than (n, TO)-dependence). 

The Mercator maps V;mooth(^, 7) = M3mooth(/3, 7; ^t) 
and Vdiscrctc (/3, 7) = V;iiscrete(a! = 0, 7; C = ;n,TO), re- 
spectively calculated via Eqs. (|2.2al) H2.2c|) and 12. 3|) 
and both based on pair potential (j2.2b|l . arc shown as 
Figs. 13 -0 for the tubes listed in Table |l| Figs. |3 Hand 
ISlare for tubes with radii around 6.5 A, 7.5 A and 8.5 A, 
respectively. Within each figure, Subfigs. (a), (b) and (c) 
refer to zig-zag, chiral and armchair tubes; the left plot 
is Vdiscrete, the right V^mooth- Only the variation is plot- 
ted; for each plot the lowest occurring energy value, for 
which we write V^", has been subtracted to make the min- 
ima lie at zero. The V'^ values for Vdiscrctc and Ksrnooth 
and the the upper bounds of the left and right plots in 
Figs. 121 - 13 exhibit discrepancies. They originate from 
the intrinsic impossibility of the smooth-tube approxi- 
mation to correctly account for the actual distribution 
of the carbon atoms on the cylinder. The wider the 
tube, the more atoms (the higher A^t), and the smaller 
the discrepancy: the (9,9) tube (V^^o^th = 70111.5 K, 
^discrete = -79117.2 K) exhibits the largest difference 
Kilcrctc - Kmooth whilc for the (21,1) tube the values 
Kfiscrctc and get very close (V^^^^^,^ = -31028.6 

K, Ki°i,„„tc = -31033.9 K). We point out that any con- 
tinuum approach suffers from such a discrepancy, and 
that it can not be resolved by replacing a in Eq. H2.3|l by 
the adjusted density a = AT/(27ri?T2zniax), a notion we 
illustrate in Appendix IXI Nevertheless, it is not senseless 
at all to perform a smooth-tube approach, because con- 
clusions are to be drawn based on the potential energy 
variation: of interest are the locations of energy min- 
ima, corresponding to molecular orientations which are 
most stable. 

The Vsmooth plots in Figs.|2Ia), (b) and (c) are, apart 
from different energy ranges, similar. They exhibit 12 
minima (white) and 20 maxima (dark gray). (The 7- 
coordinate, ranging from to 27r, is cyclic, i.e. the molec- 
ular orientations at the 7 = "edge" are repeated at the 
7 = 27r "edge". The /?-coordinate is not cyclic: points 
along the (3 = and (3 = tt "edges" with equal 7 refer 
to distinct configurations.) The 12 angle pairs 
corresponding to minimal-energy configurations are tab- 
ulated in Table Hll and indicated schematically in Fig. El 
Each of the 12 minima corresponds to a molecular orien- 
tation where two opposing pentagons of the Ceo molecule 
are perpendicular to the tube axis [Fig.j^b)]; the Euler 
transformations (a = 0, 7^), i = 1, . . . , 12, yield four 
truly different types of orientations (Fig. Hand Table ITl|l . 
The 20 maxima correspond to situations where two fac- 
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TABLE L Characteristics of selected (n, m) tubes. Tubes of all types (zig-zag, chiral and armchair) with radii Rt as close to 
6.5 A, 7.5 A and 8.5 A as possible were chosen. The angle Aq is the rotational period of the tube when performing a rotation 
about the z-axis. An expression for the calculation of the translational periodicity Az is given in Ref. 3. 
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TABLE II: Rotation angles f3i and 7i, i = 1, . . . , 12, the cor- 
responding Euler transformations 5H~^(a — 0, l3i,^i) of which 
make two facing pentagons of the Ceo molecule lie perpendic- 
ular to the 2:-axis. For small tubes (see Fig. ^ , the minima of 
Kimooth(/3, 7) occur at these 12 molecular orientations. Four 
distinct molecular orientations, labelled I, II, III and IV, are 
obtained, distinguishable by the orientation of the top {z > 0) 
pentagon (Fig. [7|l. 
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ing hexagons of the Ceo molecule are perpendicular to 
the z-axis [Fig.|2Ic)]. 

Comparing V'discrctc and 14mooth of Fig.Ofa) — (16,0) 
tube — , we see that the locations of minima and max- 
ima hardly (or even do not) differ. For the (9,9) arm- 
chair tube, Fig. 12c), the energy ranges coincide, and the 
minima locations of Vdiscrctc and T4mooth are, if not co- 
inciding, almost equal. Deviations are observed in Fig. 
ISIb) for the chiral (14,4) tube: the minima locations of 
^discrete clearly deviate somewhat from those of Kimooth; 
some even "split into two" . The deviations are small, 
however: one may write the true Vdiscrctc minima loca- 
tions as (/?■ = Pi + Af3i, 7- — ■^i + Aji). We estimate max- 
imal deviation values at \Af3i\ w 10° and |A7i| « 12°. 

For Rt « 7.5 A (Fig. the Vdiscrctc plots become 



extremely similar to the respective Vsmooth plots. For all 
three investigated tubes — (19,0), (16,5) and (11,11) 
— the locations of minima (and maxima) can be con- 
cluded to coincide. The 30 minima correspond with two 
opposing double bonds being perpendicular to the tube 
axis [Fig. Ufa)] . Maximal energy occurs when two oppos- 
ing pentagons are perpendicular to the tube axis (the 12 
minimal-energy configurations for Rt ~ 6.5 A, Fig.|31). 

In Fig. El i?T ~ 8.5 A, Vdiscrctc and V^mooth match 
completely. With respect to Fig.|21 minima and maxima 
have been flipped: lowest-energy configurations now fea- 
ture hexagons perpendicular to the tube axis, pentagons 
perpendicular to the tube axis yield the highest energy. 

Up to now, a and C have been kept fixed. To get 
an idea of the energy variation when a and C are al- 
lowed to vary, we have calculated Vdiscrctc (a, /? = 0,7 = 
0; C = 0; n, m) and Vdiscrete(a = 0, /3 = 0, 7 = 0; C; m), 
for < a < Aa and < C < AC, respectively. 
The tube-dependent rotational and translational peri- 
ods Aa and AC are given in Table HJ In Table IIIII 
we summarise these calculations by listing the differ- 
ences Aa Vdiscrctc = max(ydiscrctc(a, /? = 0, 7 = 0; C = 
0;n, m)) - min(Vdiscrete(Q;, /3 = 0,7 = 0;C = 0;n,m)) 

and A^Fdiscrotc = max (Vdiscrctc (a = 0,/3 = 0,7 = 

0; C; n, m)) - min(ydiscrete(a = 0, /3 = 0, 7 = 0; C; m)) 
conveying the energy variation. Clearly, the tubes with 
i?T « 6.5 A display fairly large energy fluctuations upon 
varying a and/or C, as one might intuitively guess from 
the jagged contours in Fig. (left). For intermediate 
(i?T ~ 7.5) and large tube radii (i?T ~ 8.5) the energy 
variations are small. 

We conclude that the smooth-tube approach works 
well for tube radii i?T ~ 7.5 A and higher, as seen from 
the Mercator maps in Figs. 0] and |S1 and the energy vari- 
ations of Table IIIII For small-radius tubes (i?T ^5 7 
A) , a systematic investigation addressing the variation of 
Vdiscrete as a fuuctiou of a, 13, 7 and Q is required. In the 
following section we perform such a study for the (16, 0), 
(14, 4) and (9, 9) and three more low-radius tubes. 
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FIG. 3: Mercator maps Vdiscroto(/3, 7) (left) and Smooth 7) (right), units K: (a) (n,m) = (16,0), (b) (n,m) = (14,4), (c) 
(n, m) — (9,9). The absolute minima values have been subtracted. 



IV. LOW-RADIUS TUBES: FULL ENERGY 
VARIATION 

Before proceeding to the full energy variation calcula- 
tion required for peapods with i?T ~ 6.5 A, we would 
like to reflect on actual small values of peapod radii. 
To our knowledge, both today's experimental and theo- 
retical situation do not show unanimity. Theoretically, 
different lower limits for a peapod's radius have been 
suggested, based on the outcome of the reaction energy 
AE in the reaction (n, m) -I- Cgo — > Geo@{n, m) — AE: 



exo- (AE < 0) or endothermic {AE > 0). Okada et al. 
|25ll2a | concluded from density- functional theory calcula- 
tions that for C6o@(", n) with 10 < n < 13 the reaction is 
exothermic, and, by extrapolating the results of 71 = 8, 9 
and 10, obtained a minimal tube radius of i?™'" « 6.4 A 
[2^ . Rochefort |23| , performing molecular mechanics cal- 
culations, set the lower limit at i?™'" « 5.9 A. From the 
experimental side, while it is still impossible to manufac- 
ture nanotubes — let alone peapods — with a given pair 
of indices (n,m), peapod samples with a narrow radial 
dispersion around a mean value Rt and good flUing rates 





FIG. 4: Mercator maps Vdiscrcto(/3, 7) (left) and T4mooth(/9, 7) (right), units K, for the (n, m) = (16, 5) tube. The minimal values 
have been subtracted. The maps for the cases (n, m) = (19, 0) and (11, 11) look very similar to the (16, 5) maps. 
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FIG. 5: Mercator maps Vdiscrcto(/3, 7) (left) and \4mooth(/9, 7) (right), units K, for the (n,m) = (21, 1) tube. The minimal values 
have been subtracted. The maps for the cases (n, m) = (21,0) and (12, 12) are very similar to the (21, 1) maps. 
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FIG. 6: Indication of the locations (/9i,7i), i — 1, . . . , 12, 
tabulated in Table Ell For tubes with _Rt ~ 6.5 A (Fig. OJ, 
these angle pairs correspond to the lowest-energy molecular 
orientations (Fig. . The contours are reproductions of the 
Vsmooth ~ 200 K contours of Fig. |2Ia). 



TABLE III: Energy variations AaKjiscrcto and A^Vdiscrcte 
(units K) when varying a or and keeping other variables 
fixed. 



(n, m) 


AcVdiscrctc (K) 


AjVdiscrctc (K) 


(16,0) 


51.3 


1079.3 


(14,4) 


643.2 


639.6 


(9,9) 


665.9 


8.3 


(19,0) 


< 0.1 


3.7 


(16, 5) 


0.4 


0.6 


(11,11) 


< 0.1 


< 0.1 


(21,0) 


< 0.1 


0.1 


(21,1) 


< 0.1 


0.3 


(12,12) 


< 0.1 


< 0.1 



(typically, 75%) can be produced at present. We mention 
a few (recent) experiments on peapods. Cambedouzou et 
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FIG. 7: When applying the 12 Euler transformations 
dl-\a = 0,A,70, with (/3i,7i), i = 1, . . . , 12, from Table 
im only four different molecular orientations the projections 
onto the {x, j/)-pIane of which are shown here, are obtained. 
Dashed (fragments of) lines have z < 0. The labeling I, 11, 
III and IV correlates with Table ITTl 



al. used a sample with Rt ~ 6.8 A. Maniwa et al. 
fitted x-ray diflFraction data on CgoQSWCNT peapods to 
simulations, resulting in a mean radius Rt ~ 6.76 A. 
Kataura et al. [3C| reported measurements on a sample 
having a diameter range of 6.25 A < i?T < 7.35 A, from 
which we calculate Rt ~ 6.8 A. The electron diffraction 
studies of Hirahara et al. 'Slj were performed on pea- 
pod samples with Rt « 7.15 A (SWCNTs from a same 
batch were used to synthesize not only Ceo-, but also 
C70 - and Cgo peapods). Kataura et al. oZ| reported 
high-yield fuUerene encapsulation, controlling the tubes 
to be "larger than the (10,10) tube", i.e. > 6.86 

A. Pfeiffer et al. [S^ inferred from Raman spectroscopy 
their three samples to have Rt ~ 7 A, Rt ~ 6.52 A and 
Rt ~ 6.505 A. From all these values one may conclude 
that peapods with a radius around 6.5 A, our represen- 
tative value for the "pentagonal case" (Figs.|21and[7Jl, al- 
though possible, are less abundant than peapods with a 
radius around, say, 6.75 A. We have therefore considered 
a few additional tubes — (17,0), (14,5) and (10,10) — 
with radii around 6.75 A; their characteristics are shown 
m Table These tubes can be expected to be more re- 
alistic representatives of the "pentagonal" regime instead 
of the Rt = 6.5 A tubes of Refs. M andlH — we recall 
that the transition from the "pentagonal" to the "double- 
bond" lowest-energy orientation occurs around 7 A [l6j |. 
We note that the (10, 10) tube is of special interest since 



tubes with a radius close to i?T(10, 10) = 6.86 A are fa- 
vorable for Cgo encapsulation, as seen in both experiment 

— according to Kataura et al. j^S] , peapod samples tend 
to have a radial dispersion centered around i?T(10, 10) 

— and theory — of the CQo@{n, n) series {n — 8, 13), 
the n = 10 peap od stands out as the most "exothermic" 
(see above) [25[l26|. 



The Mercator maps Kiniooth(/3, 7) and Vdiscroto(^, 7) of 
the extra tubes are shown in Fig. |S1 and the energy vari- 
ations AaVdiscrctc and A(;Vdiscrete are listed in Table HVl 
Again, minimal energies occur around (/3i,7i), and max- 
ima correspond to hexagons being perpendicular to the 
tube axis. The (17, 0) tube's Vdiscioto 50 K contour devi- 
ates from its smooth-tube 50 K contour [Fig. |SJa)], but 
there is some over-all agreement. The (10, 10) tube's 
Vdiscrctc plot [Fig. |Sfc), left] features "split" minima as 
seen for the (14,4) tube [Fig. Elb)]. The (14,5) tube's 
Vdiscrcto plot does not coincide nicely with its Kimooth 
plot, but interestingly, the two locations (/3i,7i), i = 3 
and 1 = 9, type III "pentagonal" orientations, correspond 
very well to their smooth-approximation counterparts. 
Since all other ten minima locations are related to the 
i = 3 or the i = 9 location by a molecular rotation over 
a = J or a = ^ about the z-axis — see Fig.[71 — , any of 
the (/3i,7i) points can be made a minimal configuration 
by changing a. The main conclusion is that the minimal- 
energy orientation will always feature facing pentagons 
perpendicular to the z-axis. The same can be said of any 
of the tubes of Figs. El and 0— tubes with i?T < 7 A — 
excepting the (14,4) tube. 



We now turn to the a- and C-dependencies of Vdiscrote- 
It is sufficient to consider the intervals < a < 
Aa = amax and < C < AC = Cmax- We di- 
vide the interval [0, Q;max[x [0, Cmax[ into a 10 x 10 grid 

(a»=i,...,io,0=i,...,io), a, = (i- l)^f^, Cj = (j- 1)%^, 
and construct a double Fourier series — for notational 
simplicity we drop the indices (n, m) — , 



Vdiscroto(a,/5,7;C) 



00 00 



Apq(P,l)cos(p^—2Tr) cos(g-^2^ 

p=0 q=0 I ^ Q^max / \ C,max 

+ Bpq{(3,^) sin[p— ^27r) cos[(7— ^27r 

\ Q^max / \ vnax 



a 



C 



+ C'pq 7) cos ( p- 27r ) sin( q- 27r 

^ Smax 

+ Dp,(/3,7)sin(p— ^27r)sin(g-^2^) L (4.1) 



an 
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TABLE IV; Characteristics of additional tubes (Rt ~ 6.75 A). See captions to Tables III and IIIII 



(n, m) 


chirality 


Rt (A) 


Aa 


AC (A) 


AaVdiscrete (K) 


AijVdiscrcto (K) 


(17,0) 


zig-zag 


6.7370 


27r/17 


VSa = 4.3128 


< 0.1 


176.4 


(14, 5) 


chiral 


6.7603 


27r 


24.5237 


36.3 


36.5 


(10,10) 


armchair 


6.8640 


27r/10 


a = 2.49 


23.1 


11.6 



^discrete (K) Kmooth (K) 




U ry- 'In Q '->|- 2n 



FIG. 8: Mercator maps ydiscrete(/3, 7) (left) and V;mooth(/3, 7) (right), units K: (a) (n, m) = (17,0), (b) (n, m) = (14,5), (c) 
(n, m) = (10, 10). The minimal values have been subtracted. 



10 



by numerically calculating the Fourier coefficients 



pq 







Cn 



da / dCVkiscreto(a,/3,7;C) 







Q^max / \ Siinax 



via the trapezium-rule approximations 

10 10 . 
. — . . — . n — r 

-Vdiscrctc(ai,/3,7;Cj) 



10 10 

i=i j=i 



C^max Cn 

^i) 10 



X cos( p 2tt I cos I q 2tt 



The coefficients Bpq, Cpq and Dpq are obtained by replac- 
ing cos() cos() by sin()cos(), cos() sin() and sin()sin(), 
respectively. The prefactors read 



^00 = 



1 



Q^maxCmax 



Q^maxCn 



Q^maxCn: 



(4.2) 



where r and s ^ 0. In series H4.1|l . some terms may vanish 
because of symmetry reasons. 

The Fourier coefficients Apg(/3, 7), Bpq{(3,j), Cpq{f3,j), 
Dpq{l3,j) can be interpreted as Mercator maps. The 
magnitude of the coefficients decreases for increasing in- 
dices p and q; we approximate Vdiscrete(aj /?, 7; C) by 



4 4 r 

Vdiscroto(a,/3,7;C) = ^P'3(/^'7) 
p=0 q=0 { 



X cos( p ^ 2it \ cosi q—^—2'n: 

Q^max / \ Simax 



(4.3) 



For given [3 and 7, we scan the a- and (^-intervals and 
define 

KiTscretel/?: 7) = min{ Vdiscrctc(a, (3, 7; C); 

0<a<amaxO<C<Cmax}. (4.4) 

The quantity ^Xscrctc(/^' 7)' gives the lowest attainable 
energy Vdiscrete when varying a and In Fig- IHl it has 
been plotted for every of the six tubes investigated. The 
plots again exhibit icosahedral symmetry as in the pre- 
vious Mercator maps. The main observation here is that 
for all tubes, except the (9,9) tube [Fig. (SJc), left], the 



absolute minima lie not precisely at the 12 (/3i,7i) loca- 
tions but somewhat away from them — the same effect 
observed for the (a = 0, C = 0) Vdiscrete plots in Figs. 
|31 and |S1 As before, we can write the actual minimum 
locations as (/?■ = (ii -\- A(3i,jl =7^-1- A'fi). For the 
Rt ~ 6.5 A tubes (Fig.El left), excepting the (9, 9) tube, 
the {(3i,'^i) locations (orientations) have energies ^ 300 
K higher than the minimal energies. The (9, 9) tube's 
minima are really close to — if not, coinciding with — 
the {Pi, 7i) orientations. For the Rt « 6.75 A tubes (Fig. 
El right), excepting the (10, 10) tube, the {l3i,ji) orienta- 
tions have energies ^ 20 K higher than the minimal ener- 
gies. The (10, 10) tube's absolute minima lie also off the 
(A,7i) "pentagon" orientations (not visible on the plot), 
and have energies ~ 12 K higher than the lowest energies. 
We must therefore conclude that chirality-dependent ef- 
fects manifest themselves here. However, for Rt « 6.75 
A tubes, probably the smallest peapod tubes as discussed 
above, the effects can be said to be minor. As an approxi- 
mation, one may consider the smooth-tube approach. We 
recall that for higher tube radii, the smooth-tube approx- 
imation is excellent (when the Ceo molecules lie on the 
tube axis). 

To conclude this section, we come back to the 4 types 
of "pentagonal" orientations depicted in Fig. [3 We recall 
that they arise from the Euler transformations D\^^{a = 

0, /3i,7i) with the angle pairs i = 1, . . . , 12, of 
Table Hll Clearly, all 12 14mooth(/3i, 7i) values are iden- 
tical because of the cylindrical symmetry. A priori, 
Vdiscrctc(ft, 7i) can be different for each of the four types 

1, II, III and IV. Interestingly, depending on the tube's 
symmetry, some orientations sometimes are equivalent. 
This is illustrated in Table Ivl and can be understood by 
being aware of certain symmetry elements. Orientations 

I and II are related via a rotation over tt about the z-axis, 
and likewise for orientations III and IV. Orientations III 
and I are related by a rotation over 7r/2 about the 0-axis. 
The presence of a twofold symmetry axis for the (n, m) 
tube therefore implies equivalence of orientations I and 

II and of orientations III and IV, while a fourfold axis 
implies the equivalence of all orientations. The (16, 0) 
tube is an example of the latter, a (14,4) tube provides 
an example of the former. The (16, 5) tube has no pure 
rotational symmetry axis and exhibits therefore distin- 
guishable Vdiscrete valucs. We note that any occurring 
energy differences as a result of the discussed inequiva- 
lences are small, however, and not seen on any of the Mer- 
cator maps since the contour values lie not sufficiently 
close to each other. Also note that these observations 
are generally true for the values Vdiscrete (A, 7i) — they 
do not need to be to minima, e.g. the (14,4) tube [Fig. 
|3fb)]. The equivalence relationships between the twelve 
(/3i, 7i) orientations is nicely seen in the the (14, 5) tube's 
^discrete plot [Fig. Elb)] : all 4 types I ({1,5,7,11}), II 
({2,4,8,10}), III ({3,9}) and IV ({6,12}) "behave" dif- 
ferently. When making the observations concerning the 
equality /inequality of 14mooth(/3i, 7i) values summarized 
schematically in Tabled we have made sure that suffi- 



11 




(16,0) Kw,. (K) 






cient numerical accuracy has been achieved. One needs 
accurate enough atomic coordinates of the Cgo molecule, 
so that upon explicitly applying coordinate transforms 
$H~^(q!, 0, 7) corresponding to symmetry elements of the 
Cfjo molecule (i.e., of the icosahedral group Ih), the same 
set of coordinates is obtained up to the desired accuracy. 
We use coordinates with 12 significant digits. 



V. DISCUSSION AND CONCLUSIONS 

We have presented a systematic comparison of the 
potential energy of a Ceo molecule — treated as an 
icosahedral cluster of ICs — encapsulated centrally in 
a SWCNT, when approximating the tube as a uniform 
cylinder and when taking the true carbon atomic net- 
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TABLE V: Schematical presentation of Vdiscrotc(/3i, 7i) values for selected tubes, indicating equalities and inequalities due 
to the presence/absence of certain symmetries. The angles fii and 7^ are tabulated in Table the corresponding molecular 
orientations depicted in Fig. ITTl 



I (type) 




Vdiscrete(/3i, 7i) 






(n,m) = (16,0) 


(n,m) = (14,4) 


(n, m) — (16, 5) 


1,5,7,11 (I) 






^a6,5) 


2,4,8, 10 (II) 
3,9 (III) 
6, 12 (IV) 


^^6,0) 


^^4,4) 
^^4,4) 
1/(14,4) 


^a6,5) 

^^6,5) 
1/(16,5) 



work into account. The former approach results in only 
two variables (the Eulcr angles [3 and 7), while the lat- 
ter requires in addition a third Euler angle a and the 
molecule's 2-coordinate C, denoting its position along the 
tube axis. The (/3, 7)-dependence can be conveniently 
plot as a Mercator map. Fixing (a,C) at (0,0) then al- 
lows a first visual comparison of the "smooth" and "dis- 
crete" Mercator maps Vsmooth (/9, 7) and Vdiscrctc(/3, 7). 
From these preliminary comparisons one can see that the 
larger the tube, the better the smooth-tube approxima- 
tion. Indeed, the Kimooth and Vdiscroto Mercator maps 
for tubes with radii i?T > 7.5 A (Figs. E| and Fig. [SJ 
are as good as identical. For smaller tubes, the effect 
of the tube structure comes into play and deviations be- 
tween the smooth and discrete Mercator maps are visi- 
ble (Figs. Eland Fig.|SI). While the (a, C)-dependence of 
^discrete can be argued to be negligible for the larger tubes 
{Rt ^ 7 A) because of the similarity of V^inooth(/3, 7) and 
Vdiscretc(/3, 7), a fuU investigation of the variables a and 
C, is in order for smaller tubes. We have presented a 
detailed study for selected zig-zag, chiral and armchair 
tubes with radii around Rt ~ 6.5 A and 6.75 A, in- 
cluding the (10, 10) tube, nowadays considered the ideal 
peapod tube [25ll26ll30| . 

A double Fourier series captures the (a, ^-dependence 
in a manageable way. Scanning Vdiscrete(a, /9, 7; C; 'ti) 
for its lowest attainable values when varying a and C, 
yields Mercator maps similar to Vsmooth(/5, 7; w). We 
see that the actual energy minima do not correspond to 
the 12 "pentagonal" (A, 7i) orientations but that they he 
slightly away from them (except for the (9, 9) tube, where 
the minimal-energy molecular configurations are really 
close to the "pentagonal" orientations). Such (/3, 7) ori- 
entations correspond to "tilted" molecules, where an axis 
connecting the midpoints of two opposing pentagons does 
not coincide with the tube's long axis (z-axis) but makes 
a cone with a small opening angle if one would perform 
the a Euler rotation. Hence, we conclude that for these 
smaller tubes (again excepting the (9,9) tube), the chi- 
rality of the tube does play a role. However, as seen from 
the small deviations occurring in the Mercator plots, the 
effect is not strong, especially for the i?T ~ 6.75 A tubes 
(Fig. 1^1 right). Therefore, one may consider the smooth- 
tube approach if one wants to capture radius-dependent 
properties. In this respect, one should be aware of the 



present-day experimental situation: precise knowledge of 
the components (i.e., occurring chiralities) in a peapod 
sample is absent — only a determination of the tube ra- 
dius distribution seems feasible up to now. The main con- 
clusion of Refs. and 0| stands: three regimes can be 
distinguished (i?T ^ 7 A, "pentagonal" orientation, 7 A 
< i?T < 7.9 A, "double-bond" orientation and i?T > 7.9 
A, "hexagonal" orientation). The "pentagonal" orienta- 
tions have to be restated as "tilted pentagonal" orienta- 
tions, though. Generally, our findings are in accordance 
with those of Troche et al. 0] who concluded that the 
chirality of the SWCNT encapsulating the Ceo molecules 
has only a minor effect. 

Note that transversal motion of the Cgo molecule (off- 
axis displacements) is not discussed here since our sole 
purpose was a direct comparison of the smooth-tube and 
the discrete-tube approaches for a centrally located Ceo 
molecule. As expected intuitively and demonstrated in 
Refs.lHandlll from a certain radius (i?T ~ 7 A) on, an 
off-center position is energetically more favorable. This, 
however, involves energy differences much larger than 
those seen upon varying (3, 7, a and C- 

The smooth-tube approximation's requiring only the 
two Euler angles (3 and 7 allows for the use of symmetry- 
adapted rotator functions. For details we refer to Refs. 
IT5I and0. The advantage lies in the possibility of writing 
Kimooth(/3, 7; ^t) as an expansion into functions il/(/3, 7) 
taking the icosahedral molecular symmetry of the Ceo 
molecule and the cylindrical symmetry of its site into 
account, 

00 

^smooth(/3,7;^T) = ^U./(i?T)il,(/3,7), (5-1) 

where the coefficients wi relate to the icosahedral sym- 
metry and carry the details of the pair interaction po- 
tential. The molecular symmetry implies the first non- 
vanishing terms to be those with I = I — I — IQ 
and 1 — 12, and a restriction to only these four leading 
terms already approximates Kimooth(/3, 7; ^t) extremely 
well. Apart from providing mathematical/physical clar- 
ity, expansion 1)5. l|l greatly reduces the calculation time. 
While, for a 100 x 100 {f3,j) grid, an implementation of 
Eq. H2.3|l takes hours, the calculation of T4mooth(/3, 7; Rt) 
via Eq. (|5.1|) is a matter of minutes on the same machine. 
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It is interesting to note that double-walled carbon 
nanotubes (DWCNTs) allow the encapsulation of Ceo 
molecules in inner spaces smaller than observed for SWC- 
NTs: Khlobystov et al. [T^ reported the insertion of Cgo 
molecules in DWCNTs with internal radii as small as 5.5 
A. Having different minimal internal radii of SWCNTs 
and DWCNTs for filling with Ceo molecules is attributed 
to the difference in how a Ceo molecule interacts with a 
SWCNT and a DWCNT |ig. We have carried out cal- 
culations for Ceo@DWCNT by treating the field on the 
Ceo molecule as a superposition of the two fields from 
the tubes with different radii. Although we find that the 
presence of a second (outer) tube decreases the energy 
for encapsulation when taking a tube radius difference 
equalling the interlayer distance of graphite (3.35 A), the 
effect is rather small and not sufhcient to explain the 
large reduction in inner tube radius. 

We believe that the general conclusion — the smooth- 
tube approximation being justified for intermediate and 
large tube radii {Rt > 7 A) and possibly acceptable 
for smaller tube-radii — reached here is relevant for 
other peapod systems. For example, (C7o)Ar@SWCNT 
peapods feature different orientations of the encapsu- 
lated C70 molecules for different tube radii |^ , the 
so-called "lying" (for smaller tube radii) and "standing" 
(larger tube radii) orientations. A smooth-tube approach 
would make a good start for investigating these specific 
orientations. 
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APPENDIX A: DISCREPANCY TEST 

We consider a single carbon atom in a short (16, 0) 
tube fragment defined by |z^| < a/\/3 w 1.42 A, leaving 
only three "rings" of 16 carbon atoms each. The tube 
atoms have coordinates rv = (x,-, yr, -^t)) t = 1, . . . ,48. 
The single atom, put at the center of the fragment which 
we define to be the origin of the employed cartesian co- 
ordinates system, has a "discrete" energy 



approximated by the "smooth" energy 



discrete 



48 



T=l 



(Al) 



V^mooth = CtRt 



2tt 



dZv^'ip). 



(A2) 



In both equations, the pair potential v^{r) of Sec. ^ 
is understood. In Eq. (|A2p . $ and Z are defined via 

TABLE VI: Vsniooth values, units K, for various tube frag- 
ments and surface densities a and a. 



[Zn, in 5 



(K) 



\/3a \/3a i 

' 2 ' 2 J 

^/3a %/3a 1 

" 2 ' 2 J 

_ 5a a 1 

~ Q 1 

4^3' V3> 
__ a a 1 

a a 1 

~ V3' 2^3^ 



a 
a 
cr 
cr 
a 
a 



-405.0 
-303.7 
-317.0 
-317.0 
-217.1 
-325.7 



X = i?T cos y = i?T sin <& and z — Z . The integra- 
tion boundaries Z^i^ and are not well-defined. In- 
deed, there is a range of both lower and upper bound- 
aries corresponding to a "smooth" tube fragment con- 
taining only the three "rings" with Zr = — a/\/3, z^ — Q 



-\/3a/2 < Zniin < — a/\/3 and 



and Zr = a/(2\/3): 
a/(2V3) < Zmax < \/3a/2. In Table EH we present 
14mooth values calculated for a few of these [Z^-m, ^max] 
intervals. For each case, both the tube surface den- 
sity a = 4/(v^a^) and the adjusted tube density a — 
48/(27ri?T(^max — ^min)) has been Considered. The "dis- 
crete" value, obtained via Eq. (|A1|1 . reads Vdiscrctc — 
—318.1 K, and is best reproduced by the "smooth" value 



if the interval [^min, Z„ 



= [-^, ^] is chosen. 



The 



tube fragment edges then lie precisely in the middle of 
two neighboring "rings" of atoms; the surface densities 
a and a then happen to coincide. The values of table 
[^suggest that some choice(s) of intervals may yield the 
^discrete value, but the point we want to make here is 
that making use of the adjusted surface density a does 
not "convert" the I4mooth to the Vdiscrctc value. We re- 
mark that doing the Vsmootii calculations described in the 
paper with a instead of a turned out to yield only very 
small differences. 
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